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Abstract 

The competition between the long range Coulomb attractive and the short range 
repulsive interaction in ionic nanocrystals creates an effective negative pressure, 
which causes the lattice parameter a to increase with decreasing nanoparticle size 
d. A simple "Madelung model" is used to predict the dependence of the lattice 
parameter for CeC-2 and BaTi03 on d, 5a/ a = a/d, for a =0.22 A and 0.18 A 
respectively. The model predictions are compared with experimental results. 
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How properties of ionic crystals change when their surface to volume ratio be- 
comes large is an active area of research [1-7]. The ceria nanoparticles Ce02 
have been studied experimentally using Raman, X-ray diffraction, photoab- 
sorption, high resolution transmission electron microscopy (HRTEM) [1,2], 
and transmission electron microscopy [3-5]. The fascinating property of ionic 
nanoparticles is the lattice parameter increase with nanoparticle decreasing 
size [1,4,5]. Intuitively one would expect the lattice parameter to decrease as 
a result of surface stress as was observed in gold nanoparticles [8]. 

Nanocrystal properties are significantly different from their bulk counterparts. 
Previous theoretical treatments [6,7] used realistic potentials for small Ce n 02 n 
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cluster calculations. It was shown that flourite structure for ceria was not 
achieved for n smaller than 50. In the present work we investigate lattice 
parameter dependence on a cluster size for d = 4 — 20 nm or n = 10 3 — 10 5 . 
Charge neutral octahedra nanocrystals are surface terminated with oxygen 
anion layers along the eight {111} faces. The balance between the long range 
Coulomb and the short range atomic interactions leads to the optimal lattice 
parameter ao dependent on the size of the cluster. The volume energy per 
formula unit in the bulk material is usually modeled as 



where B is the bulk modulus and E v = Ne v is the volume energy of N formula 
units of Ce0 2 , each filling volume space V = ajj/4. In the infinite crystal Eq. 
(1) is the total free energy which is minimum at the bulk lattice parameter a . 
In the nanocrystal the balance between the Coulomb and short range forces 
is altered due to the surface. The surface Coulomb energy creates a negative 
pressure which leads to the lattice parameter expansion in the nanocrystal in 
accord with observation [1]. 

To calculate the Coulomb energy in ionic nanoparticles we use the same pro- 
cedure as Madelung energy calculation in bulk material [9]. In the fluorite 
(CaF 2 ) structure the Ce 4+ ions occupy a fee sublattice and O 2- ions occupy 
the tetrahedral interstitial sites forming a simple cubic sublattice of length 
cio/2. This gives a Madelung energy per formula unit ^Mad = Amike 2 /a, where 
fihuik — —46.55. The value of f3 — MziZj/r is related to the Madelung con- 
stant M = 2.52 for the CaF 2 crystal structure [10], where ionic charges are 
Zi = 2, Zj = 4 and the nearest neighbor distance is r = \/3/4 in units of a. 
To calculate the surface Madelung energy we use a finite size cluster of Ce 4+ 
and 2 ~ ions terminated by the oxygen layers along the eight (111) planes. 
The surface energy is proportional to the area of nanoparticle ~ N 2 ^ 3 . The 
diameter d = (3N/27i) 1 ^ 3 ao is related to the number of Ce atoms N such that 
the sphere of equal diameter contains an equal volume NV . The Madelung 
energy of a cluster (with a uniform lattice parameter) can be calculated for 
different diameters. The results are shown on Fig. 1(a), which are well fitted 



with two parameters /Sbuik = 46.55 and [3 snr = 4.13 nm. 

The uniform lattice parameter expansion reduces the surface Madelung energy 
creating the stress in the nanocrystal estimated with Eq. (1). The total energy 
per Ce atom on the nanoparticle is 
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Fig. 1. Closed circles are model calculations of the Madelung energy on ionic (a) 
Ce02 and (b) BaTi03 nanocrystals. The Ce02 nanocrystal was terminated by 
eighth (111) planes and BaTiC-3 by six (100) planes. The solid lines are the best fits 
using Eq. (2). 
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Here the first term is the volume contribution, which includes the Madelung 
energy of the infinite crystal, and the second term is the surface Coulomb 
energy from Eq. 2. Minimization of the total energy with respect to the lattice 
parameter for fixed N gives 
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The bulk modulus B = 230 ±10 GPA has been measured by high-pressure 
x-ray diffraction [11]. Taking a$ = 5.409 A we obtain the relation 5a/ ao = 0.22 
K/d for the lattice parameter expansion shown on Fig. 2. The experimental 
lattice parameter (also shown on Fig. 2) expands with a larger slope a = 0.37 
A. Our model predicts 60% of the measured value of a. The missing 40% of 
the lattice parameter expansion might be assigned to the formation of point 
defects. These defects are oxygen vacancies and simultaneous reduction of two 
ions of Ce 4+ to Ce 3+ for each oxygen vacancy formed, maintaining charge 
neutrality. Electron spectroscopy [12] and X-ray photoelectron spectroscopy 
[3] show the presence of a Ce 3+ component in nanocrystals. Since Ce 3+ has a 
larger radius compared to Ce 4+ there may be a lattice expansion accompanying 
the replacement of (Ce 4+ ) 2 2 ~ by Ce 3+ V^ where V is a neutral oxygen vacancy. 
An upper limit in the Ce 3+ contribution to expansion was estimated [1] by 
interpolation of published results of the variation in ceria lattice parameter 
with RE 3+ radii from doping with different rare earth oxides (RE 2 03) [13,14]. 
Recent neutron diffraction measurements on nano-scale powder of Ce02 [15] 
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Fig. 2. Filled circles are experimental lattice parameters versus inverse diameter of 
CeC>2 nanoparticle. Dashed line is the best fit for the experimental data [1]. Solid 
line is a Madelung model prediction 8a/ a = 0.22 A/d. 

suggest associated oxygen defect formation (involving oxygen vacancies and 
self-interstitials) in the fluorite structure. Such defects are not included in our 
model and may alter the surface Madelung energy estimation. 



The Madelung model is valid for any ionic nanocrystal. We apply it for the 
BaTi0 3 nanoparticles to calculate the competition between the long range 
Coulomb and the short range atomic interactions. The results are shown on 
Fig. 1(b). Ionic charges Ba 2+ , Ti 4+ and O 2- were used to calculate the surface 
energy of the nanoparticle. To preserve charge neutrality the cubic nanocrystal 
with size d was terminated by the three BaO and three TiC>2 {100} planes. 
The calculated value of /3 SUT = 3.6 nm along with the experimental lattice 
parameter a = 3.996 A and bulk modulus B = 196 GPa [16] predicts 8a/ a = 
0.18 A/d. Existing data [5] on BaTi0 3 nanoparticles are not sufficient to verify 
the relationship Eq. (4) and further experiments are needed. 

In conclusion we propose a mechanism for lattice parameter expansion of 
ionic nanocrystals due to the effective negative Madelung pressure. Using the 
experimental bulk modulus we quantitatively explain a lattice parameter ex- 
pansion in ionic materials. A discrepancy between the model prediction and 
observed effect in Ce0 2 can be attributed to the Ce 3+ component [3,4] and 
other defects not included in our model. The Ce 3+ concentration deduced from 
magnetic measurements [4] suggests a much smaller value than expected from 
lattice expansion measurements in accord with our theory. A non-uniform lat- 
tice parameter expansion and structural defect formation [15] would change 
the quantitative answer, but not the conclusion that an effective negative 
Madelung pressure is significant in ionic nanocrystals. 
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